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Abstract 

A relativistic equation is deduced for the bound state of two particles, by assuming a 
proper boundary condition for the propagation of the negative-energy states. It reduces to 
the (one-body)Dirac equation in the infinite limit of one of the constituent mass. It also 
has the symmetries to assure the existence of the anti-bound-state with the same mass. 
The interaction kernel(pseudo-potential) is systematically constructed by diagonalizing the 
Hamiltonian of the background field theory, by which the retardation effects are included 
in the interaction. Its wave function at the origin(WFO) behaves properly in a manner 
suggested by the covariant field theory. The equation can well be applied in the heavy 
quark effective theory. It can also be used as a starting equation for the unperturbed 
system in the rigorous perturbation theory of QED. 



1 Introduction 

One of the unsatisfactory nature of the Bethe-Salpeter equation for the bound state 
is that it does not reduce to the Dirac equation in the infinite limit of the one of the 
constituent mass, even when the interaction is assumed to be instantaneous 0. We 
have to sum up all the crossed diagrams to recover the Dirac equation ||, which is 

*A preliminary report was presented in the workshop "Fundamental Problems in the Elementary 
Particle Theory" held at Ninon University in March 1995. 
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impossible for the finite masses. The reducibility to the Dirac equation is essential 
in applications to the heavy flavored quarkonia or the other atom-like systems. 

Historically, the relativistic single-time equation for the two-body system pre- 
ceded the BS equation. Soon after the discovery of the Dirac equation, Breit pro- 
posed the equation of the form|| 

{H 1 (p 1 ) + H 2 (p 2 ) + V}^ = E^ (1) 
where Hi is the Dirac Hamiltonian 

Hiipi) = on -Pi + rriiPi 

and V is a local potential. The Breit equation reduces to the Dirac equation in 
the limit mentioned above but does not have the ".E-parity symmetry", by which 
we mean there is symmetric negative eigenvalue E for each positive one, which is 
interpreted as the bound state of the antiparticles. 

On the other hand, in the instantaneous BS equation we have an extra projection 
factor 

A ++ - A__ (2) 
in front of the potential V, which consists of the energy-projection operators 

A £7? (Pi,P 2 ) = AKPi)A?(P 2 )> e,r] = + oi-, (3) 

where 

Al( Pi ) = {Eiipi) + eH i (p i )}/2E i (p i ), 

E t (p t ) = (p' + m^. 

.E-parity symmetry of this equation is a result of the PCT invariance: The Dirac 
Hamiltonian is odd under the PCT transformation and the interaction Hamiltonian 
is assumed to be invariant under it. We therefore see that the symmetry is assured 
by a projection factor @. We must introduce a similar factor in any attempt at the 
construction of an improved two-body equation. 

The factor (||D comes from the Stiickelberg-Feynman boundary condition for the 
propagation of the negative-energy state Q. We will construct the equation for the 
unequal-mass constituents by imposing similar boundary condition. Before pre- 
senting it, we restrict the framework of consideration. For definiteness, we assume 
Abelian gauge field interacting with the Dirac particles. We also work in the rest 
frame of the bound system(P = (E,0)), since we are looking for the non-covariant 
approximation of the low-energy dynamics, p and x, in the following, are the relative 
momentum and coordinate, respectively, in this frame. 
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2 The equation for unequal-mass constituents 



In this section, we deduce a new single-time equation (Two-body Dirac equation) for 
the unequal-mass constituents, for which we assume that the mass M(= mi) is larger 
than m(= 7712). 



2.1 Two-body Dirac equation 

Let us first consider the 2-body propagator for which we impose the boundary condi- 
tion that the negative-energy states propagate backward in time. If we keep individ- 
uality of the constituents and use the usual Feynman propagator in the instantaneous 
BS equation, the equation (P modified by the factor @ results. However, we can 
choose the other possibility in which we incorporate the idea that the bound two 
bodies should be treated as a quantum-mechanical unity.Q Since M is larger than m, 
the free part of the Hamiltonian has the same sign as that of the particle 1 in the CM 
system. Let us then modify the boundary condition as follows: A bound two-particle 
state propagates backward in time if their net energy is negative. By assuming it, we 
have the propagator 



sP(p, P ) -i: XiE _ po _ Hii _ p) - — s X2E +po _ H2{p) - — 5 ^ihi (4) 

where Ai + A2 = 1 and the limit 5 — ► +0 is assumed. 

By assuming (Q) we obtain the single-time equation. The pseudo-4-dimensional 
form of it, in the momentum space, is 

^(p) = iSP(P,p) J V(p,q)^(q)d 4 q/(27r)\ (5) 

where V(p, q) represents the interaction, which does not depend on the relative 
energies but is not necessarily an instanteneous local potential. After integrating out 
the redundant degree of the freedom, fl5|) becomes 

{ffi(-p) + H 2 (p) + J2eA ev V}tP = Eip, (6) 

£77 

which reduces to the Dirac equation in the infinite limit of M. It is easy to see the 
.E-parity symmetry of this equation. 



*We can establish the concept of individuality in the quantum mechanics only through observa- 
tion, which brings about a subtle point to the bound system: To detect an individual one in the 
bound constituents, we need to separate them by applying the 3rd interaction, which inevitably 
destroys the original state. So there is no a priori reason to apply the free propagator individually 
to each constituent. Note that "the free propagator" itself for the bound state is merely a convention 
of approximation. 
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If we would apply our equation to the scattering state, we shall have, from the 
time-dependent equation, the conserved probability density 



P (x,t) = il;(x,tyj2 £A e V ^(x,t), (7) 

erj 

which is in accord with the boundary condition that the negative- energy state prop- 
agates backward in time carrying the negative probability density 0. However, it 
does not necessarily provide the normalization condition for the bound state: For 
the scattering processes, the projected wave function K en ^ corresponds to the physi- 
cal state of the (free)particles with the positive or negative energy E. And the above 
interpretation of the probability current actually says that the state with the negative 
E is carrying the negative probability density. However, for the bound state with the 

positive eigenvalue E, A ip or A_ + ip is merely a negative- energy component in the 

representation in which the energy of the free particle is diagonal. 

Taking the above consideration into account, we restore the probability interpre- 
tation of the wave function and normalize it by assuming the probability density 

p(x)=if>(x,tyi>(x,t), (8) 

for the stationary state. Observables except the Hamiltonian are self-adjoint under 
this metric: 

(cj>,diP) = (O0,V). (9) 

The Hamiltonian is the operator ruling the time development of the system and its 
interaction part is modified by the factor £) £ „£A err Though it is not a self-adjoint 
operator, its eigenvalue is proved to be real if the inner product (|l"3"D below exists. 

2.2 Green's function and the vertex equation 

We define the Green's function G for (Q) by the operator equation 

{E — Hi — H 2 -J2eA ev V}G = $>A £ „ (10) 

erj eri 

and 

G{E-H X -H 2 -VY, = E eA ev (11) 

erj erj 

In the momentum representation, it can be written, by using the eigen-function Xnip) 
of (H), as 

G(P,P') = E TFT^ ^rXn(p)Xn(p'y + continuum, (12) 

n Nn\E — b n ) 

where E n is an eigenvalue and N n is a normalization factor defined by 

Nn = (Xn,J2 £A evXn)- (13) 

erj 
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When one of the constituents(labeled with 2) belongs to the same class of the antipar- 
ticle of the other, there can be an annihilation process through the weak interaction 
for which the unamputated-decay-vertex $ is given by 

$ = C 7 G, (14) 

where 7 is the lowest vertex and C is the charge-conjugation matrix of the particle 
2. 

$ satisfies the vertex equation 

tf 2 -V$>A e)7 ) = C 7 $>A e „ (15) 

and the amputated vertex is 

T = 7 + CW. (16) 

We can determine the renormalization constant Z\ for the wave function at the 
origin (WFO) from (fl5f) and ( pf) , if we need it||. 



2.3 Interaction Hamiltonian 

We have, so far, not specified the interaction Hamiltonian (quasipotential). In this 
section, we investigate it for the one-(Abelian)gauge-boson exchange in the Coulomb 
gauge as an example. For the instantaneous part of the interaction, V is obvious. 
For the remaining part, we can specify the quasipotential in a clear way from the 
background field theory. We have already shown, for the Salpeter equation, that the 
quasipotential from the one-boson exchange is given through the diagonalization of 
Fukuda, Sawada, TaketaniJT] and Okubo|Q (FSTO)||. We also apply this method 
to the present equation. 

We introduce the generalized Fock subspace of the free constituents, the bases 
of which are denoted by \e,7),p), where e and 77 are the signs of the energies of 
the particles 1 and 2 respectively.^ We then diagonalize the Hamiltonian in the 
Schrodinger picture by using the FSTO method. The second-order boson-exchange 
potential in this subspace is given by 



(2tt)* 



(e,rj,p\V(lb)\e',r)',p') = 77^3 X)aii(<% - -^QiQj) a 2j 



11 1 
X 2 [ <z 2 - {eE^p) - e'E 1 (p')}^ + - { V E 2 (p) - iE 2 {p')}^ 



2 There was some error concerning the Fock space in Rcf.Q. Namely, we employed the 
usual (positive-energy) Fock space and reinterpreted the matrix elements of the interaction Hamilto- 
nian including the negative-energy indices as the ones in this space according to the hole theory. It, 
however, brings the procedure into confusion, since we have revived the negative energy in Eq.(^). 
However, the error is only conceptual for the Salpeter equation which has only projection factors 



A ++ and A The result of Ref.Jl is correct. 
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where q = p — p' . The retardation effects are included in this equation. 



3 Inclusion of crossed diagrams and the counter 
correction 

We have modified the boundary condition in (^), which amounts to inclusion of 
some parts of the crossed diagram. To estimate this effect, we take up the 4th-order 
scattering amplitude in the Coulomb potential for the mass-shell particles. We are 

interested in the contribution of the intermediate state with the A H projection. 

Now, let us compare the two amplitudes restricted to this intermediate state: the 
'ladder diagram'(Z) with the propagator ([|) and the crossed diagram(c) with the 
usual Feynman propagator. We see that they are different by the following factors: 



(P2 + <?) 2 + m 2 - ctW . (p + q ) - pP) 



m 



(P2 + q) 2 + m 2 {A + s Jp 2 2 + m? + (p 2 + q) 2 + m 2 
for the 'ladder diagram' and 



(18) 



(P2 -q) 2 + m 2 -a^ ■ (p' 2 - q) - {3^ 



m 



(p 2 — q) 2 + m 2 (—A + \ p'2 2 + m 2 + J (p 2 — q) 2 + m 2 



(19) 



for the crossed diagram, where p 2 and p 2 are the initial and final momentum of 
the lighter particle respectively and q is the loop momentum to be integrated. The 
common element A is given by 



A = ^p 2 + M 2 — y/( Pl - q) 2 + M 2 . 

We see that // and I c have the same sign and the order of magnitude for momenta 
much smaller than M. Therefore, the term with the factor l\ can be interpreted as 
a substitute for the crossed diagram. 

The (pseudo)potential model is an effective theory in the strong coupling QED, 
for which the effective inclusion of the crossed diagrams is desirable. In contrast to 
it, we are able to get corrections of any desired accuracy by using the perturbation 
theory in the weak coupling QED. We could start with the Schrodinger equation 
for the unperturbed system for this purpose. However, more suitable choice is the 
relativistic equation with the good features. The two-body equation of the present 
paper is a candidate for it. 

If we intend to proceed to the next order approximation to the interaction V, 
we shall calculate the 4th-order crossed diagram by assuming ordinary Feynman 
propagator. We should, further, add counter correction terms to compensate the part 
of the crossed diagram already included effectively in the iteration of the lowest kernel 
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V mentioned above. The counter correction to the two-body propagator is given by 
taking the difference of the equation (|]) and the Feynman's two-body propagator. 
The counter corrections should be added also in the higher order approximations. 



4 Some properties 

4.1 WFO of the 1 S state 

We can use our equation as the basic equation of the effective theory of QCD. When it 
is applied to the system in which the pair annihilation of the constituents can occur, 
an important physical quantity is the wave function at the origin. For example, the 
decay amplitude of the pseudo-scalar Qq meson via a weak boson is proportional to 
the average WFO Tr{7 5 7 -?/;(0)}, where ip(0) is the charge conjugated (with respect 
to the particle 2(g)) WFO. We investigate, in Appendix, the asymptotic behavior 
of the momentum-space wave function by using the method given in Ref.||. We 
assume instantaneous exchange of a gauge boson[]. The average WFO thus obtained 
is finite. This result is consistent with consideration on the covariant field theory. 
We note that the average WFO becomes divergent in the limit of the one-body Dirac 
equation, for which we have the renormalization procedure [fLOf] . 

There are many "two-body Dirac equation" proposed. An interesting one from 



the point of view of the present paper is the one by Mandelzweig and Wallace ||11||. 
They intended to include the effects of the higher-order interaction (the crossed 
Feynman diagram) and got an equation which has the proper one-body limit and 
the .E-parity symmetry. An important difference from our equation is in the average 
WFO considered above. It is divergent in their equation if the transverse part of the 



gauge-boson exchange is included [12 



4.2 On application to the heavy quark effective theory 

One of the research fields in which we can utilize the two-body Dirac equation is the 
physics of the heavy flavored quarkonia. The recent trends in this field are led by 
the heavy quark effective theory (HQET). The Qq system is well described by using 
the two-body Dirac equation with some phenomenological (pseudo)potential. If we 
take the heavy quark limit (M — > oo) of the quark Q, the equation itself becomes the 
one-body Dirac equation, which is the basic equation of HQET. 

There are, however, divergent (as M — ► oo) portions in the perturbative correction 
to the WFO. For example, the matrix element of the current for the annihilation 
decay of a pseudoscalar quarkonium is given by 

3 The analysis in the Appendix cannot be applied to the retarded interaction. 
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< 0\q lul5 Q\Qq >= -C{u)tr{ lvl5 *(0)}, (20) 

where ^(0) is the WFO of the bound state and we have separated the correction C(fj) 
coming from the range of momenta fi ~ M in the integral over the loop momenta. 
C (n) is determined from the perturbative loop corrections summed up by using the 
renormalization group equation (T^j . 

C<J*) = £^) d , d= -6/(33 -2JV». (21) 

To improve the approximation, we can calculate 1/M corrections to HQET by 
using the two-body Dirac equation. 



4.3 Another approach to the heavy quark phenomenology 

An alternative way of investigating the heavy quark systems is due to direct appli- 
cation of the two-body Dirac equation. The binding interaction, in this formulation, 
consists of the Coulomb and the confining potentials. We further add the gluonic 
correction term to the potential, for which we should attach the high momentum 
cutoff [^]. The correction factor C(fi) depends on this cutoff. If we take the cutoff 
of the QCD scale, we should include all the corrections from the high momenta, 
except for the Coulomb contribution, to the vertex correction. By subtracting it in 
the equation (plj) , we get d = 10/(33 — 2Nf) for the exponent.^ 



4.4 On the equal-mass limit 

The unequal-mass equation (|6]) is well applied to the system with M 3> m. For 
M ~ m, the reason justifiing (f|) becomes obscure and we will have two different 
equations in the limit M = m. Though we concern ourselves in the case M > m, it 
is meaningful to investigate the degree of ambiguity near the equal mass limit. For 
M = m, the projection factor in front of the interaction term of (|J) includes a part 
which violates the exchange symmetry: It is shown that 

A( y ) = A + _ - A_ + (22) 
and the Heisenberg's exchange operator 

Ph = ~(1 + <Tt ■ ct 2 )(1 + Pl ■ p 2 )P M (23) 

4 We have assumed the counter correction which compensates the modification of the two-body 
propagator, since it cannot be regarded as a substitute for crossed diagram in QCD. If we do not 
include the counter correction we shall get d = 2/(33 — 2iV/)|Q. 
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anticommute, where the operator Pm exchanges the momenta(or coordinates). 
violates the symmetry since the remaining part of the Hamiltonian and Ph commute. 

For equal-mass limit, we have two equations. One is the equation (|B|) with M = m 
and another is obtained by assigning a minus sign in front of A^, which is the equal- 
mass limit of the equation with m > M. It is the conjugate equation of (§) in the 
sense that @ is converted into it by the transformation Pjj. It is easy to show that 
these equations have the common eigenvalue spectrum: If an eigenfunction \n of @ 
belongs to some eigenvalue E n , PuXn is the solution of the conjugate equation with 
the same eigenvalue. However, Xn does not have the definite P^-parity. 

For equal masses, it is reasonable to use the Salpeter equation which includes the 
projection factor (Q). The eigenvalue of this equation is different from the above E n . 
The difference is, however, small since it is of the 4th order in the symmetry-breaking 
part of the Hamiltonian. 

5 Summary and discussion 

By assuming a proper boundary condition for the two-body propagation in the 
negative-energy state, we have modified the propagator and proposed a new bound- 
state equation for the unequal-mass constituents. It has the symmetrical energy 
eigenvalues E n and — E n and reduces to the (one-body) Dirac equation in the infinite 
limit of one of the constituents masses. Secondly, we have discussed the normal- 
ization of the wave function and pointed out that the positive-definite probability 
density should be assumed. We can consistently calculate the observables of the 
bound state by assuming this normalization. 

The interaction Hamiltonian of the equation is constructed by diagonalizing the 
field theoretic Hamiltonian in the generalized Fock subspace of the two particles. 
Relativistic effects such as the retardation are taken into account in systematic way 
in the framework of the perturbation theory. 

We can use this two-body equation as the foundation of the rigorous pertur- 
bation theory of the weak-coupling QED. In this case, we should correct back the 
modification of the propagator. 

We have further investigated the WFO's of the proposed equation in some detail 
for the instantaneous interaction and shown that the average WFO in the % state 
which determines the leptonic decay rate of a pseudo-scalar meson is finite, which is 
in accord with the expectation from the field theory. 

One of the system for which we may utilize the two-body Dirac equation is the 
heavy flavored quarkonium. There are various possibilities of choosing the framework 
of the approximation according to the treatment for the high-momentum interaction. 
We have briefly discussed the correction factor to the leptonic decay width caused 
by the high momenta. The equation affords a good foundation of the heavy quark 
effective theory. 
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Finally, we discuss some feature of the spectrum of the equation (|6]). There are 
physical eigenvalues E n which reduce to M + m in the weak coupling limit and 
the corresponding negative ones which are interpreted as the bound states of the 
antiparticles 1 and 2. Besides those, there may be a series of eigenvalues which 
reduce to M — m in the weak-coupling limit and its negative counter part, which 
are unphysical. Fortunately, we can identify and discard these unphysical spectra by 
inspecting the weak-coupling limit. 
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Appendix 

We examine the asymptotic (p — > oo) behavior of the momentum-space wave function 
and show that the average WFO Tr{757oV>(0)}/V2 is finite. 5 

There are 4 partial amplitudes h er) (p) in the 1 S'o state, with which the wave 
function is expanded as 

X(P) =EE c r<(-p)v- r (p)h EV {p)(—^—y/ 2 , (24) 

where c\i% = — c_i/2 = 1/a/2 and the spinors u and v, for the particle 1 and 2 
respectively, are defined in |J. 

The average WFO for the annihilation decay through the axial-vector current is 
given by 



£{V(£i + eM)(E 2 + r]m) - er ] J(E l - eM)(E 2 - r,m)}h ev ( P )p 2 d P . (25) 



We first assume the Coulomb potential. The partial-wave equation for the 1 5 , 
state is given by 



{E - eE^p) - V E 2 (p)}h £ri (p) = -e— T,j d< i 

e'r/' 



'See Ref.flrn] and references therein, for the Salpeter equation. 
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+ {ee'A^Al, + m'A^A^^Q^h^q), (26) 
where z = (p 2 + q 2 )/2pq and Qe(z) is the Legendre's function. A l ££ , is defined by 



A\ £ , = y/(Ei(p) + emi){Ei{q) + e'rrn). 

The asymptotic behavior of the wave function is determined from the integral 
region near the infinity. We then expand the both sides of ( p6]) into the series of 
1/p and \jq. We assume the power behavior of the wave function for large p. The 
independent amplitudes are chosen to be h A (p) = J2 £ h ££ (p), h B {p) = J3 £ e/i ££ (p), 
hcip) = S e ^e-e(p), and hr>(p) = J2 £ £ f l e-e(p), which are expanded, in the high- 
momentum region, in power series of 1/p: 

n 

Integrals on the right-hand side can be done if we neglect infrared- divergent terms 
which are irrelevant to the leading asymptotic behavior. Now, we can determine the 
asymptotic indices /3x's from consistency @: We get, for Ha and hs, 

2C° A p-^ - EC° B p-^- 1 = -C° A —^— cot(^)p-^ (27) 

7r 1 - p A 2 

and 

2C B p-^ - EC\p-^ = -C% ^}~^\ ta< (28) 

where the terms of the higher power in 1/p are neglected. If we neglect the second 
term in the left-hand sides of (|2~7D , we find (3 a in the range 1 < f3 A < 2 [] and get 
Pb = Pa + 1 from (p8|). We obtain another series by neglecting the second term in 



fl28l) . For this, (3 B is found to be in the range 2 < p < (3b < 3, where the lower 



bound Po corresponds to the upper bound 4/n of a above which the index Pa from 



(j2"7D becomes complex. Pa of the second series is given by Pa = Pb + 1- 

The asymptotic amplitudes h c and /ip are determined dependently on Ha and 
Kb- We get, for the minimum indices 

Pc = mm(p A + 2,p B + 1) (29) 
Pd = Pa + I- (30) 
We see that the average WFO ( p5|) is finite, because 

Pb > 1 and p c > 2 

hold for the asymptotic amplitudes. This conclusion is valid even if the instantaneous 
exchange (transverse part) of the gauge boson is added. 

6 Sec Rcf.fll for the details. 
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